14 Fractals 



We now proceed to quantify the "strangeness" of strange attractors. There 
are two processes of interest, each associated with a measurable quantity: 

• sensitivity to initial conditions, quantified by Lyaponov exponents. 

• repetitive folding of attractors, quantified by the fractal dimension. 

Now we consider fractals, and defer Lyaponov exponents to the next lecture. 

We shall see that the fractal dimension can be associated with the effective 
number of degrees of freedom that are "excited" by the dynamics, e.g., 

• the number of independent variables; 

• the number of oscillatory modes; or 

• the number of peaks in the power spectrum 



14.1 Definition 



Consider an attractor A formed by a set of points in a p-dimensional space: 

-e- 




We contain each point within a (hyper)-cube of linear dimension e. 
Let N{e) = smallest number of cubes of size e needed to cover A. 
Then if 

N[e) = Ce~^, as e ^ 0, C = const, 
then D is called the fractal (or Hausdorf) dimension. 
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Solve for D (in the limit e ^ Qi): 

\nN{e)-\nC 



D = 



\n{l/e) 



Since InC/ ln(l/£:) — ^ as e: — ^ 0, we obtain the formal definition 

£-0 ln(l/£) 

14.2 Examples 

Suppose ^4 is a line segment of length L: 



Then the "boxes" that cover A are just line segments of length £, and it is 
obvious that 

N{e) = Le-^ ^ D = 1. 

Next suppose A is a surface or area S. Then 

N{€) = ^ D = 2. 

But we have yet to learn anything from D. 

Consider instead the Cantor set. Start with a unit line segment: 





The successively remove the middle third: 



1/3 2/3 



1/9 2/9^ _ 

etc etc 
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Note that the structure is scale-invariant: from far away, you see the middle 
1/3 missing; closer up, you see a different middle 1/3 missing. 

The effect is visually similar to that seen in the Lorenz, Henon, and Rossler 
attr actors. 

The fractal dimension of the Cantor set is easily calculated from the definition 
ofD: 



Obviously, 
Then 



Thus 




Taking s = 1/S and using the definition of D, 

In 2^ ln2 

D = lim = - — ~ 0.63 

m^oo In 3"^ In 3 



14.3 Correlation dimension i/ 



We proceed now to an alternative procedure for the calculation of the fractal 
dimension, which offers additional (physical) insight. 

Rather than calculating the fractal dimension via its definition, we calculate 
the correlation dimension v. 

We shall show that v < D. But first we define it. 



14.3.1 Definition 



Consider a set of points distributed on a plane. 
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Let N{r) = number of points located inside a circle of radius r. 
Assume the points are uniformly distributed on a curve like 




For r sufficiently small compared to the curvature of the curve, we have 

N{r) oc r 



or 

N{r) oc r^, u = 1. 

Now assume the points are uniformly distributed along a surface in two di- 
mensions: 




Now 

N{r) oc 

Next, reconsider the Cantor set: 



u = 2. 



We expect that N{r) will grow more slowly than r. 
Indeed, calculations show that u ~ 0.63 = D, just as before. 
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14.3.2 Computation 



Our implicit definition of v is clearly generalized by considering 

• an attractor in a p-dimensional space, and 

• N{r) = number of points in a p-dimensional hypersphere of radius r. 

For a time series we reconstruct a p-dimensional phase space with the 
coordinates 

x{t),x{t + T),x{t + 2r), ...x{t-\-{p- l)r) = x{t). 

Suppose there are m points on the attractor. We quantify the spatial corre- 
lation of these points by defining 

C(r) = lim — - (number of pairs i,j for which \xi — Xj\ < r) . 

m— >oo ui 

More formally, 

= i™. ;;;2 E E - 1^* - 

where 



^ 3 



H{x) = 



I x>0 
else. 



The summation is performed by centering hyperspheres on each of the m 
points. 

In practice, one embeds the signal x{t) in a phase space of dimension p, for 

p = 2,3,4,5,... 
p is called the embedding dimension. 

For each p, we calculate C{r). Then, assuming 

C(r) = 
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we plot logC vs. logr and estimating the slope v: 




logr 



Consider the example of white noise. Then x{t) is a series of uncorrelated 
random numbers, and we expect 

C{r) oc r^, p = embedding dimension. 

Graphically, one expect a series of plots like 

5 

log C(r) ^ ' 




log r 



Here 

u{p) = p, 

a consequence of the fact that white noise has as many degrees of freedom 
(i.e., independent "modes") as there are data points. 

Consider instead X{t) = periodic function, i.e., a limit cycle, with only one 
fundamental frequency. 

Then the attractor looks like 
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Provided that r is sufficiently smaller than the curvature of the limit cycle, 
we expect 

C(r) ocr\ for p = 2, 3,4,... 
Graphically, we obtain 

P= 2 3 4 5 




log r 

and therefore 

uIp) = 1, independent of p. 

We conclude that v measures something related to the "number of degrees 
of freedom" needed to parameterize an attr actor. 

Specifically, suppose a dynamical regime has n oscillatory modes. The at- 
tractor is then a torus T"^, and we expect 

C(r) oc r". 

Thus 

p < n =^ C(r) oc 

and 

p > n =^ C(r) oc r", independent of p. 

Conclusion: If, for embedding dimensions p > po, is independent of p, 
then 1/ is the number of degrees of freedom excited by the system. 

This conclusion provides for an appealing conjecture: since white noise gives 

iy{p) = p, 

v independent of p (and reasonably small) implies that the signal is deter- 
ministic, and characterizable by ~z/ variables. 
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(There are some practical limitations: 



• r must be small compared to the attr actor size. 

• r and m must be large enough for reasonable statistics.) 



14.4 Relationship of z/ to D 



(Orassberger and Procaccia, Physica 9D, 183 (1983)) 

The correlation dimension is not strictly the same as the fractal dimension, 
however it can be. We now derive their mathematical relation. 

Suppose we cover an attr actor A with N{r) hyper cubes of size r. 

If the points are uniformly distributed on A, the probability that a point falls 
into the ith hypercube is 

Pi = l/N{r). 

By definition, for an attractor containing m points, 

^ m m fin 

I j 

C{r) measures the number of pairs of points within a distance r of each other. 
In a box of size r, there are on average mpi points, all within the range r. 

Therefore, within a factor of 0(1) (i.e., ignoring box boundaries and factors of two arising 
from counting pairs twice), 

^ N{r) 

C{r) ~ -^^{mp,f 

i=l 

Nir) 

= 

i=l 

Then, using angle brackets to represent mean quantities, we have, from 
Schwartz's inequality, 

C(r) = Nir) (p?) > Nir) (pif = ^ 



Nir) 
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If the attractor has fractal dimension D, then 

N{r) oc r~^, r ^ 0. 
The definition of the correlation dimension v, on the other hand, gives 

C(r) oc r". 

Substituting these relations into both sides of the inequality, we find 

r" > 

Thus as r — ^ 0, we see that 

The equality is obtained when {pf) = {pi)'^. 

Thus V < D results from non-uniformity of points on the attractor. 
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